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SCATTERING  MECHANISMS  FOR 
SEMICONDUCTOR  TRANSPORT 
CALCULATIONS 


J.  Bade 
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Universiiy  of  IHinoit  at  Urhana-Champaign,  Illinou  6IS01 


Monte  Carlo  simnlatioiu  for  trniuport  ia  temieoadactoi*  anmetieallj 
solve  tlie  Boltsmaaa  eqnatioa  while  ofTeriag  the  phjrsicalljr  iatoitive  pietnze 
of  free  flights  sad  cattiet  scatteriags  oa  a  iniczoseo|Me  levd.  The  aim  of  a 
good  Moate  Carlo  simnlator  is  to  ^ve  the  most  physically  correct  tealisa* 
tioa  of  the  seatteriag  aad  free  flight  processes  ia  aa  efllcieat  maaaer.  This 
chapter  treats  the  seatteriag  mechanisms  aad  their  implemeatatioa  ia  de> 
tafl.  The  basic  scattering  meehaaisms  treated  in  this  diapter  frU  into  the 
categories  of  phoaoiu,  static  imparities,  and  seatteriag  dne  to  the  conlomb 
potential  of  other  carriers  (carrier-carrier  interactions).  Photon  seatteriag 
is  important  in  specific  eases  (radiative  reeombmatioa  rates  fiw  instance), 
bnt  occnxs  on  time  scales  mnch  longer  than  the  scattering  processes  men¬ 
tioned  above,  and  so  is  not  nsnally  a  fhctor  in  determiaiag  basic  transport 
parameters. 

In  almost  all  approaches  to  electron  transport  ia  crystals,  the  dectroas 
are  decoupled  from  the  ions,  imparities,  aad  from  each  other,  aad  <»»»»nUT 
Bloch  states  are  assumed  for  a  complete  basis  set  lot  the  dectroas.  Then, 
the  effect  of  each  of  these  interactions  is  rdatrodnced  as  a  perturbation  to 
the  simple,  one-dectron  states  forcing  tranritions  between  them,  ie.  scatter¬ 
ing.  In  order  to  accomplish  the  decoapling,  one  makes  use  of  the  adiabatic 
approximation  (Bom-Oppeahdmer  approximation)  which  is  the  topic  of 
the  first  section.  However,  the  Boltsmsuin  equation,  which  is  the  ba^  for 
all  Moate  Carlo  simulations,  is  a  fully  classical  equation,  oripaally  derived 
for  transport  of  gasses  whose  coiutitneats  interact  weakly.  The  geaeralisa- 
tioa  of  this  equation  to  the  semiclasrical  repme  ia  which  are  seen 

as  localised  events  in  space  and  time  but  calculated  quantum  mechanically 
is  a  simple  extension  of  the  dasrical  Boltsmann  equation.  The  steps  lead¬ 
ing  to  this  extenrion  which  are  rdevant  to  the  treatment  of  scattoiag  are 
discussed  in  the  second  section  of  this  chapter.  It  is  then  easy  to  obtain  the 
proper  form  for  the  scattering  rates.  Then,  a  short  discussioa  of  didectrie 
screening  (aa  important  consideration  for  the  calculation  of  scattering)  is 
given.  The  remainder  of  the  chapter  ia  devoted  to  the  calculation  of  spe^c 
scattering  rates  and  their  implemeatatioa  ia  Monte  Carlo  rimulations. 


1  The  Adiabatic  Approximation 

To  begin  ont  diKusrioa  of  election  tinaspoit  in  semiconducton,  we  define 
elections  to  be  in  one  of  thiee  classes:  core  eUeiront,  whicb  ate  tightly 
bound  to  the  nnclei,  wUnet  tUeirotn,  which  ate  loosely  bound  and  foim 
the  covalent  bonds  between  ions,  and  conduction  eleetront.  Valence  electron 
states  lie  beneath  the  band  gap  and  aie  completely  filled  at  sero  tempet* 
atuie.  Elections  which  occupy  the  exdted  states  above  the  band  gap  aie 
termed  conduction  electrons.  In  the  following  discussion  of  the  adiabatic 
approximation,  the  crystal  ts  decomposed  into  ions  (the  nuclei  and  the  core 
elections  treated  as  a  unit  moving  rigidly  together)  and  the  outer  eUetrons 
which  ate  either  valence  electroiu  or  condnction  electrons. 

The  efieci  of  the  lattice  (ions)  on  electron  tiaiupoit  is  typically  eval¬ 
uated  by  making  use  of  the  adiahatie  approximation  which  separates  the 
electronic  (outer  electrons)  and  ionic  motion  so  that  their  interaction  may 
be  treated  in  perturbation  theory.  The  crystal  Hamiltoniaain  the  adiabatic 
approximation  can  be  derived  from  the  general  Hamiltonian  as  follows.  The 
full  Hamiltonian  is  (1), 

BtM  =  T|«.({H})-Hf,«.({i2})-|-r.,({r})  +  H„({r})  +  H.,({r}.{H})  (1) 

where  the  ionic  and  electronic  kinetic  energies  are  Tttm.  uid  Td,  and  the 
potential  energies  Utom  and  H**.  (U,,  is  the  electron-electron  interaction.) 

is  the  electron-ion  interaction,  and  the  ion  and  electron  coordinates 
are  denoted  by  {R}  and  {r}.  Because  the  ions  are  more  massive  than  the 
electrons,  they  will  move  much  more  slowly.  The  key  step  in  formulating 
the  adiabatic  approximation  is  the  assumption  that  the  electrons  respond 
adtohatically  to  the  motion  of  the  ioiu  —  the  ion  motion  does  not  force 
tranritions  between  electronic  states,  but  the  electron  eigen-states  evolve 
adiabaticaDy  as  the  ion  positioiu  change.  Then,  for  a  fixed  set  of  ion  coor¬ 
dinates  {J{},  there  corresponds  a  frdl  set  of  electronic  eigen-functions: 

mi(W)  +  «^..(W)  (2) 

and  the  eigen-function  of  the  full  Hamiltonian  pven  by, 

♦=  E  *•«*»♦•«'}■  w)  w 

[T,^m)  +  +  E^]  Xn({R})  =  X-({il})  .  (5) 

In  (5)  several  small  terms  have  been  neglected  which  involve  the  action  of 
on  the  electronic  wave-function.  This  is  a  justifiable  approximation 
since  these  terms  have  been  shown  to  be  of  order  (m/Af)*^’  where  m  is 
the  electron  mass  and  Af  is  the  ion  mass  (see  [1]).  The  ionic  wavefnnetions 
must  be  calculated  self-conristently  with  the  dectronic  wavefnnetions  as 
shown  in  (S). 

In  order  to  untangle  the  dection  states  from  each  other,  the  Hartiee  or 


Hutree-Fock  approsmation  can  be  employed  which  decoaples  the  many* 
electron  state  into  simple,  one  electron  states  by  the  introduction  of  ex* 
change  and  correlation  enerpes.  The  electron-electron  interaction  Utt  ia- 
troduces  another  effect  which  dresses  the  remaining  coulomb  interactions, 
such  as  Ufi,  with  a  screening  cloud  resulting  in  screened  interactions 
section  3).  We  shall  designate  this  screened  electron-ion  interaction  as  U,{. 
Similarly,  the  electron-electron  polarisation  screens  other  conlomb  interac¬ 
tions  such  as  the  interaction  of  conduction  electrons  (holes)  with  ionised 
impurities,  K,  and  the  electron-electron  (carrier-carrier)  interaction  be¬ 
tween  conduction  electrons  or  between  conduction  electrons  and  valence 
electrons,  Vt,.  K  and  are  assumed  weak  and  treated  as  the  perturba¬ 
tions  responsible  for  ionised  impurity  scattering  and  impact  ionisation  (the 
inverse  Anger  process)  discussed  in  later  sections  of  this  chapter. 

The  many  electron  state  is  retrieved  by  placing  electrons  into  the  one 
electron  states  obeying  the  Pauli  exclnsion  principle.  This  is  the  essence  of 
the  one-electron  approximation  for  which, 

[Tw(r)  +  {«})  -I-  U,x]  W) .  («) 

where  Uax  i*  tlia  exchange  energy. 

The  electron  mobilities  are  generally  large  in  most  semiconductors,  and 
it  is  appropriate  to  treat  the  electron-ion  interaction  as  a  perturbation 
(2).  Furthermore,  because  ionic  vibrations  involved  are  relativdy  small,  a 
good  choice  for  a  pertnrbation  approach  to  decouple  the  electrons  fiom  the 
ions  is  an  expansion  in  the  displacements  of  the  ions  from  thcix  equilibrium 
positions  «*(iZ|).  If  these  eqnilibrinm  positions  ate  designated  as  {JZ*’},  the 
solution  of  (6)  for  the  one-dectron  wave-functions  yields  the  famihat  Blodi 
wave-functions  of  band  index  n  and  wave-vector  k: 

[T.,(r)-hI^.i(r.OT)-l-irnx]l^-»{r)  =  (7) 

^•(r)Mr)  =  Mr)  =  i«-»(r)e«-'  . 

Labeling  the  equilibrium  lattice  sites  by  JIj  for  integers  /,  the  basu  vectors 
o,  and  keeping  only  the  first  two  terms  in  the  expanrion  of  —  J{|), 

^.<(r  -  fl,)  «  &Mr  u-(Jl,)  •  VxMr  -  M  ^ 

the  one  electron  crystal  Hamiltonian  in  the  adiabatic  approximation  be¬ 
comes: 

=  T,,  +  V,(r)  +  V,(r,J{)  +  V,{r)  +  V„(r,r')  (9) 

y,(r)  = 

Here,  Vt(r)  is  the  equilibrium  periodic  crystal  potential,  &ti(r,  {A**})*  l^(r) 
represents  the  phonon  scattering  perturbation. 

The  ion  eigen-states  are,  (to  second  order  in  the  ion-ion  interaction) 


hftrmonic  o«cill«tot  sUtet  or  phoaoaa.  SpedfleaU/,  the  lea  Ramiltoaiaa 
u  diagoaalired  by  appealiag  to  the  taieiag  aad  loweiiag  operators  for  the 
aotmal  modes  of  vibratioa,  the  phoaoa  modes  of  the  crystal.  In  the  phoaoa 
oecnpatioa  aumber  basis,  {3] 

+  <(«)«/(»)).  (10) 

ti 

where  at  (g)  aad  aj(q)  ate  the  tatsiogaod  lowetiag  operators  for  the  phoaoa 
mode  of  wavereetor  q  aad  braaeh  j.  la  this  aotatioa,  caa  be  ex> 

pressed  by  aa  ezpaasioa  over  the  aotmal  modes  of  wareveetor  q  of  the 
lattice  aa  [4] 

•*(»)  =  E  + •/ (-1)1 4  (Ji) 

where  N  is  the  anmber  of  primative  cdls  ia  the  lattice,  nia  is  the  mass  of 
the  ioa  at  bads  loeatioa  a,  aad  ^  is  the  phoaoa  waTcyeetor. 

2  Scattering  in  the  Semi-classical  Boltzmann 
Equation 

To  detiye  the  ezptesdons  for  scatteriag  ia  the  seoii-clasaical  Boltsmaaa 
eqnatioa  (SCBB),  we  start  with  the  eqaatioa  of  motioa  ibc  the  deasity  map 
triz  p  for  the  adiabatic  czystaL  As  discnssed  ia  sectioa  1,  the  crystal  system 
eaa  be  separated  iato  a  carrier  part  (eoadaetioa  dectroas  ov  holes)  aad  a 
phoaoa  pari.  The  phoaoos  are  coondered  to  be  a  thermodyaaaiie  heat 
bath  ia  eqoilibriom  at  some  temperatare  T.  Thea,  weahly  coaled  to  it  are 
the  carriers  which  we  caa  caomerate  with  the  Bloch  waveTector  aad  baad 
iadex  quaatnm  aombers.  Furthermore,  the  carriers  ate  weaUy  eoapled  to 
each  other,  K«i  aad  to  imparities,  Vj.  As  aoted  above,  these  bteractioas 
are  weak,  aad  we  caa  approximate  the  full  deasity  asatriz  as  a  product  of 
the  carrier  deoaty  matrix  aad  the  equilibtium  phoaoa  deasity  matrix  at  all 
times.  TUs  derivatioa  igaores  so  that  nm^e  oae-partide  states  caa  be 
used;  however,  to  iodude  to  flrtt  order,  two  partide  aatisymmeteriied 
states  caa  be  chosea  iastead.  For  oae-partide  states,  the  caxriet  deadty 
matrix  is  labded  only  by  the  waveveetor  k  (supprcssiBg  a)  aad  the 
phoaoa  bath  deadty  matrix  by  p#  s  Here,  fi  =  where 

kg  is  Boltsmaoa’s  coastaat,  and  Bg  =  Hi,  the  bath  (ionic)  Haaultoniaa 
given  above.  la  general,  the  quantum  analogue  of  the  phase  space  distribu¬ 
tion  function  is  the  \i^gner  distribution  which  ia  a  transformation  of  p.  To 
demonstrate  the  form  of  the  scattering  kemd  (‘coDidoa  integral*)  in  the 
SCBB,  the  Add  can  be  ignored  aad  it  is  permisdble  to  ideatity  the  diag- 
ond  of  as  the  distribution  fuactioa  f{k)  which  occurs  ia  the  dasdeal 
Boltsmann  equation. 

Rrst,  the  standard  quantum  Liouville  equation  neglecting  the  dectrie 


field  U: 


(12) 


where  the  h7pet>oper«tot  noUtioa 


A*B  =  lA,B]=:AB-BA 


(13) 


has  been  employed  and  A  and  B  ate  opetators.  Fm  is  the  adiabatic  Hamil¬ 
tonian  as  derived  in  section  1. 

The  projection  operator  method  (see  for  example  Knbo  [5])  separates 
the  diagond  of  the  system  density  matrix  so  that  an  equation  for  the  time 
evolution  of  s  f(k)  is  obtained.  Then,  the  two  projection  operators 


F 

Q 


1-P 


(M) 

(16) 


with  the  properties  PP  s  P,  QQ  =  Q,  and  PQ  =  QP  =  0,  can  be  used  in 
order  to  derive  equatiotu  of  motion  fi>r  the  diagonal  and  oiF-diagonal  parts  of 
p  in  terms  of  each  other.  Here,  tra  is  the  many  body  trace  over  the  phonon 
bath  states.  The  projection  operator  P  projects  the  portion  of  p  which 
is  diagonal  in  the  carrier  quantum  numbers  (containing  the  probability  of 
occupancy  of  the  carrier  stationary  states)  and  leaves  the  phonon  bath  in 
its  equilibrium  configuration. 

By  substitution,  we  arrive  at  an  equation  for  only  the  diagonal  part  in 
the  form, 

=  -ji-W"  jC  (i«) 

where  the  density  matrix  has  been  assumed  to  be  diagonal  at  f  a  0.  Here, 
all  the  perturbations  discussed  earlier  have  been  included  in  Bp. 

Then,  noting  the  hyper-operator  identity  for  operators  A  and  B, 

€*’b  =  e-'He-'*  (17) 

and  that 

e*'-(‘-*^  =  lf.(<,0  (W) 

where  H’«(<,  t*}  is  the  Schroedinger  pcture  propagator  for  the  full  HamHto- 
nian  from  f  to  t*,  we  can  write  equation  (16)  in  terms  of  these  propagators 
(111. 

^Pp{t)  =  J*deu,{t,i‘){s,rpp(nu,iif,i).  (19) 

Since  Bp  is  oonadered  weak,  it  can  be  neglected  in  the  exponential  of  (18) 
compart  to  Ba  und  Bo  —  T,i  +  V,.  This  allows  an  expression  of  U,  in 


terms  of  products  of  the  &ee  propsfstors  for  the  both  ond  the  carriers, 


ta  (20) 

If  the  matrix  element  for  the  electron  phonon  interaction  giren  in  (10) 
is  defined  as, 

(V,)„.S<«V|K,(r)|«h>,  (21) 

where  the  phonon  operators  in  a*(i2{)  have  produced  to  Njq  as  discussed 
above,  and  the  matrix  element  for  the  non-phonon  scattering  terms  is  de¬ 
fined  as 

(V,)„.S<nV|V|(r)|«h>  (22) 

for  the  perturbation,  then  the  equation  of  motion  for  s  /(Jfc)  is  given 

by, 

^/(*.<)*  .  (M) 


I  I’  -3i(h,t»;<0/(*.<)K  + 

A  Jn 

with 

ri(4,*')  s  (24) 

3>(*,*0  =*  +  l)e-‘<‘’t*)'*'<*'>'*'‘-*'X‘-0  + 

,  (25) 

with  w(h)  =  ^(h)/lt.  T<(h,h')  and  Tfik^V)  are  functions  associated  with 
non-phonon  and  phonon  scattering  from  kiok*  respectivdy.  Furthermore, 
the  trace  over  the  phonons  and  the  action  of  the  rairing  and  lowering  opera¬ 
tors  in  the  electron-phonon  perturbation  have  produced  the  factors  Njf  -f  1 
and  Nj^  for  emisaon  and  absorption  respectivdy  with  Njf  being  the  Bote- 
Eirutein  average  occupation  number  for  the  phonon  of  wavevector  q  and 
phonon  branch  j: 

If  /(<)  were  outride  of  the  time  integrals,  the  integral  over  t  of  the 
T(h,  V)  Hornes  the  familiar  delta  functions  of  energy  conservation.  For 
example, 

lim  f*di'Re(Ti{k,k*;f))^icKS{E{k)-E(k')).  (27) 

Since  /(t)  typically  changes  very  little  in  the  time  it  takes  to  establish 
a  delta  fuction,  /(t)  can  indeed  be  taken  outride  of  the  time  integral. 
The  result  is  the  standard  Boltsmann  equation  with  the  scattering  rate 


S(k,  k*)  bom  k  to  V  given  by  tbe  fntniiiu  GoMea  Role  ezpteaaoa  fot  energy 
conservation: 


Siik,  4')  =  j-l  (V0„.  \*S{B{k)  -  E{k'))  (28) 

for  elastic  scattering  processes,  and 

for  pbonon  scattering  processes  wbere  tbe  upper  sign  is  taben  for  emisrion 
and  the  lower  sign  fot  absorption  and  g  s  4  —4'.  Tbas,  the  scattering  cates 
in  the  SCBE  ate  rimply  given  by  the  Fermi  Golden  Rnle  bom  qnantnm 
mechanics. 

However,  the  scattering  events  have  been  treated  as  if  they  happened 
instaotaneonsly  as  a  resnlt  of  the  approamation  in  (27).  This  a^coziina- 
tion  breaks  down  in  three  eases.  The  Best  case  is  obvions.  If  transieats  ate 
to  be  resolved  on  a  time  scale  of  the  order  of  the  time  it  tabes  to  estab¬ 
lish  (27)  and  the  scattering  rates  ate  high,  /(4,t)  cannot  be  palled  oat  ci 
the  time  integral  and  the  simple  forms  in  (28)  or  (29)  will  not  apply.  If 
the  pertnrbations  ate  sufficiently  weak,  then  /(4,t)  can  still  he  removed 
bom  the  Ume  integral,  but  (28  and  (29)  mnst  be  replaced  with  partially 
completed  delta  functioiu,  which  relax  energy  conservation.  This  is  jnst  a 
manifestation  of  the  energy  and  time  ancertainty  principle. 

Second,  if  scattering  rates  ate  high,  the  ^>ptorimatioa  of  a  weak  inter¬ 
action  lea^g  to  equation  (18)  breaks  down,  and  in  the  rime  it  takes  to 
establish  the  delta  function  in  (27),  the  orig^al  state  can  decay  ^pteciably. 
To  account  for  the  depletion  of  the  initial  state  in  the  ftnite  rinm  it  takes 
fot  the  ctdlirion  to  become  "complete”  the  approodmarion  in  (18)  must  be 
improved  to  include  the  elTeet  of  in  17,(1, f').  One  way  to  account  for 
Hp  is  tc  calculate  the  full  propagator  I7«(t,f^  bom  idd  theory.  This  leads 
to  two  additions  to  (28),  First,  the  simple  propagators  in  the  scattering 
terms  ate  dressed  by  vbtual  transitioas  in  the  adf^ergy  S(4,H)  [6].  The 
second  addition  unfortnnatdy  leads  to  the  inclnsioa  of  terms  wh^  do  not 
tesemUe  the  Boltimann  equation.  Typically  these  terms  ace  assumed  small 
and  ace  ignored.  At  present,  it  is  undear  what  effect  these  terms  have.  Fbl- 
lowing  the  general  approach  that  these  terms  are  small,  we  can  augment 
the  SCBB  by  induding  £(4,  S)  in  equation  (18). 

Assumbg  that  the  self-energy  can  be  calculated,  (see  lor  instance  refer¬ 
ence  (7]  for  a  sdf-consisteat  calculation  of  the  sdf-eaergy  tot  teaSstic  band 
structures)  quantum  fidd  theory  states  that  {oti>f  [6], 

17,(4;  1,0  =  €-*<*W+*(*)K«-Oe-*(n‘)X*-0,  (jo) 


with 


A(4)  *  Re(E(4)) 
r(4)  =  /m(E(4))  . 


(81) 


Phjriieallj,  A(A)  cottespoaJa  to  a  tklA  ia  tk«  caotgj  krd  ^(A),  aad  r(A) 
eottetpoadf  to  tli«  ftaito  lifetime  of  tbe  etate  ft.  Pot  example,  the  ptobabilit/ 
that  a  pattiele  wbieb  etattc  at  time  t*  ia  etate  ft  ie  ia  etate  ft  at  time  t  (I  >  t*) 

(32) 

The  total  eeatterias  rate  oat  of  etate  ft,  5i«t(ft),  caa  be  ideatilled  ae 

SiMik)  =  |r(ft),  (M) 

bp  iatetpretiag  the  lifetime  of  the  etate  ae  the  iavetee  of  the  total  eeatteiiag 
rate.  (Tbie  caa  alto  be  ebowa  tbtoagb  applicatioa  of  the  Optical  Theorem 
qaaatam  mcebaaies  (8].)  Tlue  eleaxtp  demoaetratee  the  eoaaeetioa  betweea 
high  eeatteriag  rates  aad  the  fiaite  lifetime  of  the  etate. 

-  Wbea  (32)  le  eabstitnted  iato  (23),  the  limit  b  eqaatioa  (27)  becomes, 

A(E,£0a  (34) 

Jo 

_ TiE).+  T{^) _ 

“  ^  («  +  A(«)  -  E*  -  A(E»))^  +  (r(S)  +  r(^))*  * 

Therefore,  oae  wap  to  add  the  effects  of  high  scatteriag  rates  to  the  SOBS 
is  to  replace  S{S  —  £*)  with  A(B,  which  is  the  fkmiliar  LorcatsiaB  Eae 
shape.  The  effect  of  broadeaiag  the  ddta  fiuetioa  throa^  Ugh  scatteriag 
rate  is  termed  celKsioa  ftroodentaf. 

Lastlp,  if  a  large  electric  8dd  is  preeeat,  the  carriers  caa  be  acceletated 
appredaUp  doiiag  the  time  of  the  coOisioa.  This  is  nUtntd  to  the  «a<fw- 
eoUmouol  fitU  tffeet  la  geaeral,  this  effect  caa  also  broadea  the  eaecgp 
eoaserriag  ddta  fhactioa,  howeTer,  it  Is  more  dUBcalt  to  deal  with  ia  a 
compact  wap.  The  iaterested  reader  is  referred  to  the  folkwiag  tefereaees 
for  a  detailed  treatmeat  of  this  effect  aad  cdlisioa  broadeaiag:  [9]>[11]. 
Also,  for  aa  orerview  of  attempts  to  iaclade  coOisioa  broadeaiag  aad  the 
iatracdOeioBal  field  effect  ia  Moate  Carlo  Simalatioas  see: 

3  Dielectric  Screening 

As  discassed  ia  seetioa  1,  the  scatteriag  poteatials  Vt,,  Vt  aad  1^  are 
edf*coarieteat,  tcreeaed  poteatials  wUch  resalt  from  bare  periarbatioas. 
Becaase  ia  most  eases  the  bare  poteatials  are  simple  Coalomb  poteatials, 
v(q)  St  l/(4e'eo4*),  the  easiest  wap  to  calcalate  the  seIf<coasisteat  pertarba* 
tioBs  is  to  screea  the  bare  peitarbatioas  with  the  dielectric  fiuetioa  which 
connects  the  two. 

Suppose  the  potential,  l^a(r,t)  is  iatrodneed  as  a  bare,  eztcnal  potea* 
rial  to  the  crpstaL  The  resolting  (tree)  potential  felt  at  time  t*  aad  poritioa 
r*  is  Vui/ti*)  end  the  iadneed  darge  is  pi«(r,t).  If  1^.  is  weak  eaoagh 
that  first  order  pertarbatioa  theorp  is  adequate,  thca  we  caa  calculate  the 


ledistribatioa  of  cbatge,  **  *  lincu  funetional  of  the  self-coanjteat 

potential  Kr(r,t).  In  genetal  [18] 

=  y<l*f'y’*'x(nr';<-OMr'.0  (35) 

vhete  xCi*!  f'i  f  ~  0i  eleettonie  anaceptibility,  u  calenlated  from  8zst  ot> 
det  pettntbation  theory  in  whieh  Vt,  it  the  perturbation.  A  related  quantity 
eall^  the  dieleetrie  fonetion  ean  nmilarly  be  defined: 

Vi.(r.<)  =  J^rfJdftiry:t-f)Vu(r^,f).  (35) 

The  'inverse*  of  the  dielectric  function  in  the  quantity  we  want  fince  it 
expreuee  the  true  potential  in  termt  of  the  bare  (external)  potentiaL 
The  Fourier  tranaform  p(q  -(■  G,u)  of  equation  (35)  ean  be  written, 

Ptm(9  +  0,u)  »  f  d(q'  +  <?')x(f +  +(ril^)y*m(^  +  <»*i«)i  (37) 

Jvti 

where  ^  it  a  reciprocal  lattice  vector,  q  it  a  wavevector  in  the  first  BtiDouin 
tone  (BZ),  and  Vol  it  the  crystal  voturne.  Because  x  1*  lattke  translation* 
ally  invariant,  ^(^  >0  =  x(^  +  r*  -i-  and  the  btegral  in  equation 
(35)  beconcs  a  summation:  (14) 

^•(«  +  Ci«)  =  53x(«  +  0',fi  +  G'.«)Vr,(q +  <?,«).  (38) 

«* 

A  rixnilar  result  holds  for  c(f  +  G,q  +  <?',v): 

V„(q  +  <?,w)  *  j;c(q  +  <?,q  +  G'.w)  V^iq  +  <?,«)  .  (39) 

O' 

AppOeation  of  Poisson's  equati<m  connects  Yt,  srith  Ka  Fonder 
transfonninf  the  Potsson'a  equations  for  the  total  charge  p  s  pm  +p«s  *Bd 
the  external  charge  p,,  we  have, 

(q  +  <7)>K*,(4  +  G,«)  =  -lp(,  +  <7,«)  (40) 

*0 

(fi  +  G)>K.(4  +  G,w)  =  -lp.,(q +  <?.*,).  (41) 

Solving  (40)  and  (41)  for  pr.  and  substituting  this  into  (38)  identifies  the 
dieleetrie  function  in  (39)  as: 

t{q  +  G,q  +  G',v)  =  tao’  -  ^  x(9  +  9  -f  G* (^3) 

The  random  phase  approximatioa  (RPA),  a  type  of  first  order  perturba¬ 
tion  calculation  including  the  temperature  through  the  distribution  function 


/,(*),  define*  x 


X*{9  +  G,q  +  <?,u) 


V'  /■•(*  +  <)  ~* 

B»>{k  +  <)  -  Bn{k)  +  liht 


<  nk  I  «-‘i*+0)  '  I  n‘(k  +  q)  ><  n'(4  +  q)  |  |  «*  >  .  (43) 


This  together  with  (42)  expresses  the  effect*  of  screening  on  n  here,  externnl 
potential.  Then,  the  inverse  of  the  matrix  e(q  +  G,  f  +  G^,  where  the  rows 
and  columns  are  labeled  by  G  and  (7*.  expresses  the  true  pertarbations  Vtt, 
Vi  and  Vf  in  terms  of  the  bare  Conlomb  potentials  which  give  rise  to  the 
perturbations.  A  plot  of  e(f,g)  is  shown  in  figure  1  as  a  function  of  w  fee 
silicon  at  seto  temperature  calculated  using  equation  (43). 

Several  simplificatioru  of  this  result  are  applicable  fer  scattering  in  semi* 
conductors.  First,  the  off-diagonal  ternu  can  usually  be  ne^ected  fee  non* 
phonon  scattering  dnee  they  ate  usually  smaUet  than  the  diagonaL  They 
must  be  kept  fer  phonon  scattering  since  they  are  necessary  to  fulfill  certain 
sum  rules  ([2]).  Furthermore,  since  there  are  usually  many  more  electrons 
in  the  valence  bands  than  there  are  electrons  (holes)  in  the  conduction  (vn- 
lence)  bands,  /11(h)  fa  1  for  valence  bands  and  f»[k)  m  0  fiM  conduction 
bands.  This  is  the  seto  temperature  approximation. 

One  simple  way  to  include  the  effect  of  the  free  conduction  band  dee* 
trons  in  the  T  =  0  model  u  to  add  the  susceptibilities  fer  the  T  s  0  case, 
X*  and  the  susceptibility  for  a  fkee  electron  gas  x*  with  density  equal  to 
the  dendty  of  conduction  band  dectrons.  Two  expressions  fer  x*  given 
by  the  Thomas-Fermi  screening  theory  and  the  Lindhard  screening  theory 
for  a  free  dectron  gas  [17].  For  example,  in  the  Thomas*Fetmi  theory 


X*{9) 


.dMS,B,) 

dB, 


(44) 


where  /o  is  an  equilibrium  Fermi  distribution  in  the  conduction  bands,  and 
Bf  is  the  Fermi  energy.  Then,  a  natural  screening  length  can  be  defined  as 


€^dMB,Bf) 

9B,  ’ 


(45) 


so  that  ^ 

«(»)  =  («) 

Having  derived  rdations  between  the  bare  perturbing  potentials  and 
the  screened  potentials  which  result,  we  can  explicitly  calculate  important 
scattering  matrix  elements  for  use  in  the  SCBB  of  section  2. 


REAL  PART  OP  DIELECTRIC  FUNCTION 


IMAGINARY  PART  OP  DIELECTRIC  FUNCTION 


Figoxe  1:  Real  and  Ima^nary  paita  of  the  firequencj  and  warcTectoi  de¬ 
pendent  dielectric  function  in  silicon  calculated  using  the  random  phase 
apptoidmation  (RPA). 


4  Phonon  Scattering 

4.1  Phonon  Perturbation  Potential 

In  tills  section,  the  phonon  peitozbntion  is  rewritten  in  terms  of  the  bare 
potential  and  the  dielectric  function,  and  the  matrix  elements 
from  (21)  ate  calculated.  From  section  1,  the  pertozbing  potential  was 
found  to  be 

y,{r)  =  x;  ^  («) 

la 

where  Otilr  —  Jtj  —  a)  is  the  self*consistent  psendopotential  (effectire  po¬ 
tential  outside  of  the  ion  itself)  felt  by  an  electron  at  r  due  to  the  ion  at 
site  Jti  and  basis  vector  a,  and  Vjt,  *  &ti(r  —  Rt)  is  the  local  gradient  of 
that  pseudopotential  taken  with  respect  to  the  ionic  location.  Ptoceeding 
as  in  Vogl  [2],  we  Fourier  transform  Vf{r): 

»;(*•)  = 

lO 

where,  9  is  a  vector  in  the  first  Brillouin  tone  and  is  a  reciprocal  lat¬ 
tice  vector.  The  summation  above  must  contain  all  vectors  of  reciprocal 
space  because  Vf(r)  depends  on  the  displacements  tt‘(lZi)  which  are  not  in 
general  periodic  in  a  primitive  lattice  translation.  Since  the  displacements 
are  small,  we  assume  that  U  related  to  the  unscreened  ionic  potential 
•a  by  the  inverse  dielectric  function  €*‘‘(r,r')  as  described  in  section  S.  In 
this  context  the  unscreened  ionic  psendopotential,  «ai  >*  potential  of 
the  nudens  and  the  core  dectrons  assodated  with  the  bads  ion  a.  If  the 
crystal  is  monoatomie,  Pa  is  independent  of  ot.  »«  behaves  like  a  Coulomb 
potential  at  long  distances,  but  has  a  tepnldon  for  short  distances  due  to 
the  core  dectrons  (see  (17]).  A  mote  exact  expression  for  this  potentiJ  Is 
pven  in  the  chapter  by  Fischetti  and  Higman. 

In  terms  of  the  bare  pseudopotential  and  the  ^dectrie  function, 

dV{q  +  G)  =  ^ 

y  «  JvM  la 

,  f  d*r'e~*(r,  r')Vj*,»a(r'  -  Hi  -  «)  1 

Jvtt 

where  Vof  is  the  aystal  volume.  The  resulting  Fourier  components,  ^V(g-|- 
G),  of  the  screened  perturbation  are  referred  to  as  psendopotential  pertur¬ 
bations. 

From  section  3  we  can  write  equation  (50)  in  reciprocal  lattice  space  as, 

dV{q  +  G)  =  -i  53«“(lli)- 

la 

53€-*(«  +  <?.«  +  G')(«  +  <r)e*t^+*>-<*+®'>*a(fl  +  O') .  (50) 


The  result  of  equetion  (SO)  comes  from  the  recogaisktion  that  Vi(4*a(r  — 
-  o)  s  —  Jti  -  a),  and  an  application  of  the  shifting  and 

derivative  properties  of  the  Fourier  transform.  We  do  not  need  to  consider 
the  frequency  dependence  in  (50)  since  in  the  adiabatic  approximation,  the 
ions  are  assumed  to  be  stationary. 

Substituting  the  expression  in  equation  (11)  for  u"{Rt)  into  equation 
(50),  we  can  write, 

(q  +  G')e'H«  +  G,  q  +  G')e«(q  +  f?')e-*(«+®')- 

I 

Concentrating  on  the  sum  over  I,  we  can  make  several  simplificationa.  First, 
<?'•  JZ(  is  always  2r,  so  the  &  in  the  iiist  exponential  can  be  ignored.  Second, 
if  q*  ^  q+G"  with  G"  being  any  reciprocal  lattice  vector,  the  sum  vanishes 
because 

.  (62) 

I 

Hence,  we  can  considerably  simplify  equation  (51)  and  write  it  as 

(q  +  <r)€-»(q  +  <7,q  +  G')v.(q  +  <?)€*'••«,  (53) 

where  the  fact  that  Q  a  Vol/i^f,  the  unit  cell  volume,  has  been  used. 
Equation  (53)  is  substituted  into  equation  (48)  to  give  the  toll  deetron* 
phonon  perturbation  Vp{r)  in  terms  of  the  bare  ion  psendopotential  and 
the  didectric  function. 

4.2  General  Electron-Phonon  Matrix  Elements 

To  cdculate  the  matrix  dements  of  the  dectron-phonon  interaetioa  as  re¬ 
quired  in  the  scattering  kemd  (see  section  2,  equation  (29)),  we  need  only 
^eulate  the  matrix  dement  of  t*^**°^'  between  Bloch  states.  Pot  exam- 
ple, 

<  I  V,{r)  I  n*  >=  53  ^^(q  -K?)  <  nV  |  e'<*+®)-’  |  n*  >  .  (54) 

fO 

where  |  nh  >  are  the  normalised  Bloch  states  i^~l«nft(r)e^*''  of  band 
index  n  and  wave  vector  k.  The  phonon  raising  and  lowering  operators 
in  equation  (53)  have  already  act«l  on  the  phonon  occupation  states  to 
give  Nfj  as  in  section  2,  equation  (26).  The  dectronie  part  of  the  matrix 


element  provides  conservation  of  the  total  crystal  momentum  as  follows: 

<  nV  I  '  1  nk  >= 

jf  <•».(»■)«•*(>■)  (65) 

=  Sk-.k±,+a>  (56) 

Ja 

Here,  the  integral  in  (55)  was  factored  into  a  sum  of  integrals  over  the 
primative  cell.  The  result  of  (56)  follows  from  the  identity  pven  in  (52). 
The  G'  in  the  Ktonecker  deita  refers  to  a  general  reciprocal  lattice  vee* 
tor,  not  necessarily  the  particular  G  in  the  integral.  When  (T*  ;£  0,  the 
scattering  process  is  called  an  Umklapp  process  (assisted  by  a  reciprocal 
lattice  vector).  In  the  following  sections,  always  implies  that  Umk* 

lapp  processes  are  allowed  -*  although  the  G^  will  be 

suppressed  in  the  notation.  If  an  Umklapp  process  oeeuts,  then  an  extra 
term,  **,  appears  in  each  integral. 

4.3  Phenomenological  Phonon  Scattering  Processes 

The  form  of  the  electron-phonon  matrix  element  as  given  in  (54),  written 
in  terms  of  microscopic  quantities  is  exact  within  the  rigid-psendo>ion  ap¬ 
proximation  (see  the  chapter  in  this  book  by  Fisehetti  and  Higman)  [19]. 
However,  because  it  requites  a  detailed  knowledge  of  both  the  didectric  ma¬ 
trix  and  the  bloch  wavefunctiou,  it  is  very  diflicnlt  to  calculate.  With  some 
approximations,  nmple  forms  for  these  interaction  matrix  elements  which 
correspond  to  wdl  known  phenomenolopcal  scattering  processes  in  senu- 
condnctors  can  be  demonstrated  [2].  In  particular,  the  acoustic  and  optical 
deformation  potentials  and  the  ^ar  optic  interactions  can  be  extracted 
from  (64),  and  the  approximations  leading  to  these  familiar  processes  can 
be  examined. 

Ikom  equation  (53)  it  is  apparent  that  Vf{r)  can  contain  both  long 
range  and  ihoti  range  components.  Long-range  components  vary  neglipUy 
witlun  the  unit  cdl  and  produce  fidds,  whose  average  over  many  cell  lengtlm 
does  not  vanish.  The  long-range  components  are  therefore  tesponsiUe  fbr 
the  macroscoine  fidds  produced  in  the  crystal,  which  have  be«  identified 
with  the  phenomenological  pdar  optic  interaction  in  polar  materials,  and 
the  piesodectric  interaction  in  both  polar  and  nonpolar  materials  ([2]).  In 
contrast,  the  short-range  components  of  the  perturlnng  potential  involve 
rapid  osciOstions  within  the  unit  eeU,  and  consequently,  their  average  on 
any  macroscopic  length  scale  vanishes.  They  do,  however,  contribute  to 
scattering  through  the  phenomenological  deformation  potential  interaction. 

The  separation  of  Vp{r)  into  short  and  long-range  components  has  been 
shown  by  Vogl  by  considering  the  behavior  of  9V(g  +  G)  as  f  -•  0  [2].  We 


quote  only  the  lesnlt, 


■+E 


ja-r  <*‘*(g  -f  G,  i) 


0^0 


dViq) 


0^0 


(57) 

(58) 


whete  V^(i‘)j,  wwi  ^^(j*),,  me  the  long  and  short  range  components  respee* 
tively. 


4.9.1  Phenomenological  Deformation  Potential 

The  electron~phonon  matrix  element  for  the  short*range  interaction  is  ob¬ 
tained  by  combining  equations  (58),  (53)  and  (55).  The  xesnlting  matrix 
element  can  be  recast  into  the  form  of  the  deformation  potential  interac¬ 
tion,  and  win  be  shown  to  involve,  lot  small  q,  dilation  for  optical  modes 
and  elastic  strw  for  acoustic  modes.  For  large  q,  the  form  of  the  phe- 
nomenolopcal  intervalley  scattering  matrix  element  will  be 

We  start  by  putting  <  n'4'  |  l^(r)^,  |  nk  >  into  a  more  manageable 
form  for  a  crystal  with  a  basis  of  two  (a  =  1  and  a  s  3).  For  a  crystal 
with  a  basis  there  ate  two  phonon  types,  optical  and  acoustic,  and  they 
behave  quite  differently.  For  acoustic  inodes,  limy— o<>»(9)  -♦  0,  correspond¬ 
ing  to  a  ti^d  displacement  of  the  entire  lattice,  whereas,  for  optic  modes, 
limy— 0  <>*(9)  “♦  ^  0  corresponding  to  rigid  displacement  of  the  two  sub- 

lattices.  The  latter  case  can  cause  transitions  between  the  ad^batic  Bloch 
states  leading  to  scattering,  but  the  former  ease  cannot  lead  to  scattering 
[20].  For  small  9  it  is  possible  to  write  the  phonon  polarisatioa  vectors, 
^  *  TTmi^/^^y  and  ss  where  (jf  is  a  unit  vector, 

IS  a  constant,  and  the  plus  (minus)  sign  is  for  acoustic  (opti¬ 

cal)  modes.  Furthermore,  in  this  limit,  there  axe  three  acoustic  and  three 
optical  branches  corresponding  to  two  modes  polarised  transversely  to  the 
direction  of  propagation  and  one  poplarised  longitudinally. 

<  I  K,w„  I «»  >  =  (i»,.(r)  +  i  ±  i)  [-^]  *{,„ 

-i  8sM±y  <  »V  I  «-*•  ' ViV(0  I  «*  > 

with  V*  ipven  by  comparison  with  (58)  and  (53).  /i  can  be  chosen  to  be 
the  reduced  mass.  Since,  the  final  form  for  (59)  is  phenomenological,  the 
constant  can  be  chosen  to  accommodate  the  choice  of /i. 

To  compare  (59)  with  the  phenomenological  forms  for  deformation  po¬ 
tential  scattering,  we  expand  in  a  TWylor  series  b  q  and  look  for  the 
higest  order  non-vamshing  terms. 

For  small  q,  the  exponential  in  equation  (59)  can  be  expanded  as  e***'  m 
l+»9>r.  For  the  acoustic  phonos  case,  the  highest  order  term  ia  (59)  which 


does  not  vnnish  is  the  f  •  r  term,  sinee  ss  discntsed  nbove,  the  consUat 
term  coizesponds  to  a  rigid  displseement  of  the  entire  crystal  which  cannot 
scatter  [3].  For  a  particular  mode  y,  (59)  then  becomes 

<  «'*'  1  I  nk  =  ^yj,(T)  +\±\y 

<  n'*'  I  (fl.r,  +  +  g,r,)V;,(r)  |  nJ;  >  («0) 

Equation  (60)  can  be  put  in  the  form  of  the  familiar  deformation  po¬ 
tential  interaction  as  first  expressed  by  Shockley,  which  for  a  giren  phonon 
branch  u  [21]*[23], 

(«) 


The  terms  in  (60)  correspond  to  the  following  terms  &om  (61): 

Ea  =  <n'*'|((V/')<f,  +  (V/'),r,]  |n*>  (62) 

where  (l^'')i  is  the  component  otVj*,  and 

Su  =  5(^i«/+<i«)  (“) 

Since  the  displacement,  «,  is  proportional  to  1^**^*%  Sa  ia  the  strain 
tensor  (25).  Thus,  Ea  is  the  deformation  potential  tensor  whidk  couples 
the  locd  strain  set  up  by  the  acoustic  phonon  to  the  scattering  matrix 
clement  as  given  in  the  phenomenological  theories  of  scattering  by  acoustic 
phonons.  Equation  (61)  is  valid  for  smaQ  q,  which  corresponds  to  iatravaOey 
scattering  for  low  energies  and  represents  the  anisotropy  of  the  coupling 
constant.  For  larger  q,  Ugher  order  terms  ia  tg  •  r  come  into  the  int^ral 
in  (60),  and  hence  the  coupling  to  the  phonon  wavevector  involves  higher 
rank  tensors. 

Now,  returning  to  (59)  for  the  ease  of  optical  phonons,  the  highest 
order  term  which  can  1m  non-iero  b  the  scro-order  term.  In  contrast  to 
the  acoustic  mode  case  ,  the  sero-order  term  can  be  nonseto  because  when 
f  — »  0  for  an  optic  mode,  the  two  sublattices  are  rigidly  displaced  with 
respect  to  one  another,  and  this  can  scatter  carriers,  llins,  for  the  case  of 
optical  phonons,  we  look  to  highest  order  at  the  sero-order  term  and  write 

(64) 

-i  (lfd(T)  +  5*5)*  I  yjW  I  «*  >  -6* 

This  is  valid  for  small  q  (intravalley  scattering).  Equation  (64)  can  also  be 


wtitten  in  turns  of  an  optical  deformation  potential  rector  D,f  as, 

<  1  I  nJi  >«,  =  (65) 

The  form  of  the  electron-phonon  matrix  element  for  small  q  and  optical 
modes  involrcs  a  direct  dilation  of  the  local  crystal  strnctnre  as  seen  fiom 
the  dot  product  conpling  Dof-i.  Thus,  68  yields  the  phenomen<Jopcal  opti¬ 
cal  deformation  potential  matrix  element  ([25]).  It  has  been  shown  that  for 
symmetry  reasons  the  intravalley  optical  deformation  potential  scattering 
vanishes  for  X  and  Ft  valleys  for  the  suo-ordu  matrix  element  discussed 
above  [26]-[27].  In  genual,  the  highu  ordu  tensor  coupling  (acoustic-like) 
can  be  non-seto. 

While  for  small  q  the  phenomenolopcal  form  of  the  deformation  poten¬ 
tial  interaeUott  for  acoustic  and  optic  modes  is  quite  diiTuent,  for  large  q, 
the  coupling  contains  many  complicated  highu  ordu  tensor  modu.  The 
standard  treatment  for  these  situations  assnmes  that,  for  a  particular  intu- 
valley  tranation  (lor  example,  X-X,  F  -  X),  q  is  confined  to  a  small  cone  of 
allowed  directions.  Because  the  angular  orientation  is  fairly  constant  within 
this  cone,  the  tensor  coupling  can  change  only  a  small  amount  ibr  any  scat¬ 
tering  into  the  cone.  Then,  it  may  be  approximatdy  correct  to  treat  the 
coupling  for  utervalley  transitions  by  a  constant.  Usually  the  intuvalley 
deformation  potential  matrix  element  is  written  phenomenolopeally  as  [24] 

<n'»'|V,(r)„|«t>h=  + 

(•7) 

Both  acoustic  and  optical  modes  have  been  shown  to  participate  in 
intervalley  scattering  with  this  type  of  scattering  rate.  It  has  been  shown, 
for  instance,  that  to  fit  experimental  valuu  of  conductivity  in  siHcoa,  it  is 
necessary  to  include  acoustic  phonons  in  utervaDey  seattcriag  [28].  It  is 
evident  &om  the  nature  of  the  integral  in  (59)  that  the  coupling  is  widely 
determined  by  the  orientation  of  the  final  vaUeys  and  the  phonon  branch. 
Thus,  for  each  phonon  branch  and  each  set  of  initial  and  final  valleys,  a  new 
inter  valley  conpling  constant,  D^,  must  be  chosen,  ly^cally,  the  values 
for  us  selected  to  but  fit  experimental  data  for  a  pven  Monte  Carlo 
simulation  [29]-[Sl]. 

Finally,  lor  the  phenomenological  deformation  potential  interaction,  some 
authors  have  factored  out  the  overlap  integral, 

/(nV,ni)  =  /  d*r  ms(r)w;.j,(r)  ,  (68) 

Ja 

from  the  matrix  elements.  This  is  only  valid  for  the  long  range  interaction, 
as  discussed  in  the  next  section,  and  can  readily  be  seen  from  (56)  and  (57). 


4.3.3  The  Pheoomenologicel  Polar  Optical  Interaction 

In  this  section,  the  phenomenological  point  optic  phonon  mattiz  element 
will  be  derired  from  the  long-range  matrix  element  in  (57),  using  the  equa¬ 
tions  of  macroscopic  electrostatics  [2].  The  first  step  is  to  factor  out  of  (57) 
all  terms  proportional  to  the  Coulomb  potential  v(q)  ec  l/g*  so  that  the 
behavior  of  the  mattiz  element  for  small  q  is  clear.  These  terms  can  be 
factored  out  pving. 


<  I  I  n*  >r* 

K,(«) 


-ir«(g)/(n'*',ni) 

>  fa 


(69) 


with, 

««(r)H(»r„(T)+i±i)‘[j^]  .  (70) 

The  quantity  Z%  is  the  remainder  of  (57)  after  the  singular  l/q  terms 
are  removed  (see  VogI  [2]).  The  whole  matrix  element  is  proportional  to  1/g 
and  so,  acts  quite  differently  than  the  deformation  potential  matriz.eleffient. 
Here,  e(q,  g)  =  e(g),  the  macroscopic  dielectric  function.  Equation  (69) 
is  a  fairly  intractable  formula,  since  the  function  H*is  generally  difficult  to 
calculate.  However,  their  exists  a  simple  expression  for  Z%  which  can  be 
seen  by  unng  macroscopic  electrostatics.  This  is  only  possiUe  because  of 
the  long-range  nature  of  the  perturbation  discussed  earlier. 

The  term  containing  the  product  can  be  treated  as  the 

effective  phonon  polarisation  vector,  «.  In  adffition,  ^^(g)  i*  the  true 
macroscopic  potential  perturbation  set  up  by  this  phonon.  As  such,  it  can 
be  examined  using  macroscopic  electrostatics  and  related  to  the  polarisatioa 
wave  of  the  phonon. 

Two  equations  from  macroscopic  electrostatics  directly  apidy:  first. 


V  .J?(r)  =  V.(eo^(r)  +  i»(r))  =  0  <=>  E(g)  =  -P(g)/eo  (71) 

in  the  absence  of  excess  charge  (charge  neutrality  is  not  necessary  for  the 
proof  but  it  is  a  convenient  assumption),  and  P(g)  is  the  polarization. 
Second, 


V.V(r)  =  -eE(r)  <=►  iqV(q)  =  -eEiq)-i  (72) 


where  V(t)  is  the  potential  energy  felt  by  an  electron.  Together  these 
equations  give 

kVm(q)  =  ti-P(q)l€t)  (73) 

which  can  be  substituted  into  (3.25)  to  give 


<  n'k'  I  Vf{r)f^  I  nk  — 


I{n'V,nk)4-eP{q) 

*oq 


(74) 


Comparing  74  with  09  yields  the  polarisation  P(q)  in  terms  of  microseopie 
quantities. 

'  For  small  q,  P{q)  is  the  dipole  polarisation  set  up  by  the  optical  phonon 
mode  in  the  crystal.  For  larger  q  P(q)  acquires  higher  order  poles  (quadtupole 
polarisation,  etc.  ).  We  are  only  interested  in  the  small  q  limit  for  polar  op¬ 
tic  phonon  scattering,  so  for  these  purposes,  P(q)  ^  Pdi^att-  Fortunately, 
the  dipole  polarisation  associated  with  a  long  wavelength  mode  is  easily 
calculated  fiom  a  self-consistent,  0rst  order  lattice  dynamics  theory.  The 
result  is  that  (32], 

4  (75) 

where  e^  is  the  effective  dipole  charge  associated  with  a  longitudinal  optical 
phonon  and  u  is  the  ionic  displacement.  Furthermore,  an  application  of  the 
well  known  Lyddane-Sachs-Teller  relation  relates  the  effective  charge,  e^, 
to  the  experimentally  known  quantities  €„  «oo  the  lonptudinal  optical 
phonon  frequency  at  sone  center  uio  [52]: 

(.1)>  =  P‘) 


This  allows  us  to  rewrite  (69)  for  small  q  in  a  more  transparent  form 
involving  experimentally  known  quantities: 


<nV|V,(r),Jni>^= 


/(n'k'.nJt)  (77) 


or,  in  the  famiUar  Ftohlich  interaction  form, 

l<  I  V,ir)„  I  nk  s  l’(n't'.nJfe) 


Pure  transverse  modes  do  not  scatter  because  for  them,  4  *  ~  0* 

So,  only  longitudinal  modes  participate  in  polar  optical  phonon  scatteriag. 
AIm,  although  not  stated  explicitly,  the  Fkohlich  interaction  (75)  vanishes 
in  non-polar  crystals.  The  reason  for  this  is  obvious  on  physical  gtounds- 
in  a  non-polar  crystal,  no  dipole  polarisation  can  arise  from  optical  mode 
displacements.  Mathematically, 


^  ^  .  (h  =  -(h  (7®) 


so  the  sum  over  a  vanishes. 

The  lowest  order  (in  q)  long  range  acoustic  mode  process  which  doesn’t 
vanish  is  the  scattering  due  to  quadrupole  polarisation.  Piesoelectrie  scat¬ 
tering  results  from  quadrupole  polarisation  and  can  be  present  in  both  polar 
and  non-polar  crystals  [25]. 


4.3.3  Phonon  Scnttering  R«te« 

In  the  first  put  of  this  section  the  general  electron-phonon  matrix  element, 
and  the  standud  phenomenological  electron-phonon  matrix  elements  have 
been  presented.  To  make  nse  of  these  rates  in  a  Monte  Culo  simulation, 
we  need  to  calculate  fiom  these  matrix  elements,  the  total  scattering  rates 
out  of  a  puticnlu  state,  and  the  difierential  scattering  probabilitjr  between 
a  given  final  and  initial  state  for  um  in  the  SCBE. 

In  section  2,  it  was  stated  that  for  weak  perturbations  and  low  electron 
energies,  the  Fermi  Golden  Rule  (28)  adequatel/  described  the  scattering 
rate  into  a  puticnlu  set  of  final  states  from  a  given  inirial  state.  The 
general  form  for  the  total  scattering  rate,  5ia<(nJk),  is  the  Gcdden  Rule  rate 
summed  over  all  final  states: 

5io.(«*)  =  E  T  d(E..(Jb)-E.(k)TH)  •  (80) 


where  is  the  phonon  associated  with  the  puticnlu  scatterer  in  Vf{r). 
This  expression  can  be  calculated  exactly  using  the  full  electron-phonon  ma¬ 
trix  clement  in  (53)  as  discussed  in  the  chapter  by  Fischetti  and  Higman. 
However,  it  is  standud  practice  for  Monte  Culo  simnlations  to  take  the  - 
much  rimpler  route  and  calculate  the  scattering  rates  for  the  phenomeno¬ 
logical  matrix  elements  of  the  previous  sections. 

In  fact,  in  order  to  obtain  closed  form  expressions  for  these  scattering 
rates,  it  is  necessuy  to  assume  rimple  analytic  forms  for  the  band  structure, 
the  most  general  of  which  is  the  non-puabolic,  ellipsoidal  band  structure 
describing  the  bands  neu  the  minimum  of  a  *valley*.  The  analytic  form 
most  often  used  is, 

a  . 

where  kr  is  the  minimum  of -the  valley.  Accordingly,  this  description  is  only 
valid  neu  the  minimum  of  a  puticnlu  valley.  In  general,  the  band  struc¬ 
ture  is  very  complicated  [33],  and  in  the  case  of  hot  electron  transport,  the 
electrons  ue  Ihx  these  minima.  One  way  to  deal  with  this  problem  is  to 
calculate  the  matrix  elements  discussed  above  for  a  full  band  structue.  An¬ 
other,  much  simpler  way  is  to  calculate  the  various  scattering  mechanisms 
valid  neu  the  mininu  and  uormalise  the  scattering  rates  to  the  density  of 
states  for  higher  energies  (30].  This  is  a  reasonaUe  approach  rince  most 
scattering  rates  ue  proportional  to  the  density  of  states. 

In  this  section  the  standud  scattering  rates  valid  neu  valley  minima 
ue  discussed.  Since  the  derivatioiu  leading  to  the  scattering  rates  from  the 
phenomenological  matrix  elements  ue  well  known  and  given  in  many  texts 
(see  for  instance  [34]  and  [36]),  we  quote  only  the  results  here.  Topical 
values  for  the  phenomenolo^cal  puameters  and  material  puameters  for 
silicon  ue  pven  in  the  appendix  of  tlus  chaptu. 

However,  u  a  prelimeninuy,  we  list  bets  useful  for  deriving  them.  First, 


to  traasfoim  •  lummation  over  V  in  (80)  into  »  mote  manageable  integral 
over  k'  we  use  the  transfotmation  [34] 

<“) 

Note  that  ipta  is  conserved  in  phonon  collisions.  Furthermore,  if  the  rm 
ate  equal,  these  integrals  are  normal]/  done  in  spherical  coordinates.  Then 
the  magnitude  of  k  can  be  transformed  into  an  integral  over  final  energ/ 
Ef  by  use  of  the  transformation. 

If  the  masses  are  unequal,  the  Herring* Vogt  transformation  is  useful  because 
it  maps  the  problem  into  a  “starred-space”  which  has  a  spherical  E(k) 
relation  [36]: 

4;— (84) 

The  problem  is  then  solved  in  the  *starred-spaee*  and  then  transformed 
back  to-  the  physical  space.  As  far  as  total  scattering  rates  are  concerned, 
this  only  has  the  effect  of  replacing  m*  in  the  result  by  mjf  S  (maWiynij)^^*. 
For  intra-valley  optical  deformation  potential  scattering: 

*  8^^  K-  + 1  i]  p(H  ±  A«^) ,  (86) 

where  f(E)  is  the  density  of  states  pven  by, 

F(^)  *  (1  +  2aH)[H(l  +  oH)]»/>  .  H>0  (86) 

and  is  the  optical  deformation  potential,  p  is  the  crystal  density  and 
is  the  relevant  optical  phonon  energy.  To  determine  the  final  state  after 
scattering,  we  only  have  to  enforce  enagy  conservation,  E^k')  =  ^(4) 

±  fivgf.  Any  state  on  this  energy  conserving  surface  is  equally  probable. 
For  inter- valley  deformation  potential  scattering: 

+  5  ±  I]  ±  -  Ar,) .  (87) 

where  W(,  is  the  intervalley  phonon  energy  for  a  particular  phonon  branch, 
and  Eg  is  the  number  of  equivalent  final  valleys.  For  instance,  in  Silicon, 
there  are  six  X-mimima,  and  two  “different”  types  of  intervalley  scattering. 
For  scattering  across  to  the  X-mimimnm  on  the  same  axis  (g-scattcring). 
Eg  s  1;  for  scattering  to  one  of  the  minima  on  the  plane  perpendicular  to 
the  initial  state  axis  (f-scattering).  Eg  = 

A|«  is  the  energy  difference  between  the  initial  and  finial  state  min¬ 
ima.  For  X-X  scattering  in  silicon,  =  0.  For  F-L  scattering  in  GaAs, 
ns  0.2eV.  As  in  the  intra- valley  optical  deformation  potential  scattering 
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+  E'/^)lKaT 
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Absorption  s  0 

*m-.  =  AE:^'\b:'‘*  +  E^t^VKaT 
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rat«,  g{E)  is  the  density  of  states  in  the  final  valley.  The  density  of  states 
plays  an  important  role  in  all  scattering  rate  determination,  but  it  is  most 
apparent  in  the  expressions  tot  the  optical  and  intervalley  deformation  po¬ 
tential  cases.  The  final  state  is  selects  randomly  from  all  states  conserving 
energy,  and  the  final  valley  is  picked  randomly  &om  the  equivalent  valleys. 
A  list  of  typical  intervalley  phonons  for  Si  is  given  in  the  appendix  to  this 
chapter. 

For  intra-valley  acoustic  deformation  potential  scattering:  Be¬ 
cause  of  the  complicated  tensor  form  for  acoustic  conpling,  there  is  an  in- 
trinsic  anisotropy  in  the  coupling.  However,  as  shown  by  ConwdQ,  the  effect 
of  this  anisotropy  is  small  [24]  and  can  be  removed  by  a  snitable  averaging. 
This  amounts  to  defining  an  average  sound  velocity  ««  s  (2«2 + v!)/S  where 
and  are  the  velocity  of  sound  propagated  by  transverse  a^  lonptn- 
dinal  modes  respectively  and  an  effective  conpUng  strength  S.  B  is  usually 
chosen  to  fit  experimental  data.  The  aconstie  phonon  scattering  rate  is[M]: 


/  “*’«fajV±(,)(i  +  2offT2aJrnr»)a* 


(W) 


and  the  dimensionless  limits  of  integration:  where  El  S  . 

For  polar  optical  scattering:  The  total  scattering  rate  Ibr  polar  optical 
phonon  scattering  has  been  calculated  by  Boardman-Faweett-Swain  which 
takes  into  account  the  overlap  integrals  in  (78)  [37].  These  integrals  can  be 
evaluated  for  k  and  k*  near  the  bottom  of  non-parabolic  bands  using  k  •  p 
perturbation  theory.  The  result  u: 

I^(k\k)  =  (aftas>-f  ctct'cosfi)’  (90) 

o»  =  [(l-|-off(*))/(l-|.2off(h)]*'* 

c»  =  [off(h)/(l-|.2aff(*))]‘'* 


(91) 


where  $  u  the  angle  between  k  and  V.  Ueing  this  ezptessioa  for  the  orerlap 
integral,  the  total  scattering  rate  for  polar  optical  phonon  scattering  can 
be  written: 


where 


A  ss  {4(1 +  ««)(!  + «£')(! +  2<»^)(l  +  2a^))-* 

B  =  {2(l  +  a£)(l  +  a^)+a{7(^)  +  7(£)l)> 

£•  s  B-huto 

Becanse  of  the  1/g’  dependence  of  the  matrix  element  and  the  angular 
dependence  of  the  overlap  integrals,  polar  opticol  phonon  scattering  is  a 
anisotropic  scattering  mechanism.  To  choose  the  final  state  we  can  use  the 
Ton  Nenmann  rejection  method  to  pick  the  an^e  i  between  k  and  k*,  which 
is  pven  by  the  probability  density, 


P(a\M  -  (ata>»  4-  Oc>»  cosfi)*  nnfi  d» 

'  ^  “  yCr')  +  yiB)  -  27>/>(^)7i/*(^eos#  * 


(93) 


The  asimuthal  angle  is  completely  random  since  the  scattering  j^ba^ty 
density  is  independent  of  Thus,  we  can  ehooee  ^  with  a  nniform  random 
nnmb«  f  :  r  €  (0, 1]  by  ^  ~  2aT.  The  magmtude  of  the  final  state 
waTereetor,  k',  is  selected  by  «a«g7  eonserration  for  the  ^vea  scatteriag 
event  (phonon  emission,  or  absorption),  and  thus,  the  final  state  vector  is 
comfdetely  determined. 


5  Impact  Ionization 

The  multiplication  of  carriers  by  impact  ioaisation  is  of  central  importance 
in  the  theory  of  semiconductor  devices  both  as  a  limitiag  mechanism  and  as 
a  basu  of  device  functionality.  Impact  ionisation  b  a  tim  dectroa  process, 
corresponding  to  the  exact  inverse  of  the  Anger  process:  a  highly  energetic 
conduction  band  dectron  collides  with  a  valence  band  dectroa  which  b 
ionised  over  the  band  gap,  leaving  two  conduction  dectrons  and  a  hole.  The 
process  can  also  occur  for  holes,  in  which  a  highly  energetic  hde  creates 
two  holes  and  an  dectron.  Impact  ionisation  for  hdes  can  be  seen  as  a 
"mirror  image*  of  the  impact  ionbatioa  for  dectrons,  so  we  will  treat  only 
the  dectroa  Impact  ionisation. 

As  in  section  1,  the  dectroa>dectron  interaction  b  designated  as  V„(r,  r*), 
and  foQowiag  sectfoa  3,  we  screen  the  base  dectron^ectson  interaetioa  with 
the  didectrie  screening  function.  We  wiD  neglect  the  off-diagoad  ternu  of 


+  G,q  +  G';w)  utd  allow  q  to  take  on  ail  values  in  the  teeiprocal 
lattice. 

The  bare  electron*«lectron  intetaction  is  of  coune  just  the  coulomb  in¬ 
teraction  between  two  electrons.  Figure  2  shows  the  impact  ionisation  pro* 
cess  schematically.  States  three  and  four  are  final  conduction  band  electron 
states.  State  one  is  the  initial  electron  state,  and  state  two  is  the  initial 
valence  band  electron  state.  The  crystal  momenta  and  band  indices  are 
designated  as  hi  and  nr,  with  t  €  [ltd].  Figure  2  shows  the  electron  in 
states  1  and  2  interacting  via  the  screened  Coulomb  potential  to  generate 
electrons  in  states  3  and  4.  The  final  hole  state  corresponds  to  the  missing 
valence  electron  in  state  two. 


Figure  2:  Schematic  Representation  of  the  screened  electron-electron  in¬ 
teraction  corresponding  to  impact  ionisation.  Time  is  plotted  vertically. 
Notice  that  the  interacUon  is  retarded  due  do  dynamic  screening  effects. 

Since  Vt,{r,  rQ  is  a  two  particle  interaction  explicitly  carrying  two  spa¬ 
tial  coordinates,  its  matrix  elements  are  between  two-particle  states.  Since 
the  electrons  are  fermions,  the  two-particle  wave-functions  mast  be  anti¬ 
symmetric  linear  combinations  of  two-particle  states  with  the  Pauli  spin 
matrices.  We  write  these  two-pactide  states  as, 

-  (94) 

!^»,»,(ra)l>„»,(ri)<ri(sj)<ra(si)]  , 

where  the  w  are  the  Pauli  spin  matrices,  and  sj  and  sa  are  the  spin  co¬ 
ordinates.  The  coordinate  wave-fnnetions  ^  are  normalised  Bloch  wave- 
functions.  Also,  we  write  the  subscript  A  to  signify  the  anti-symmetrised 


state. 

la  this  language,  the  event  shown  in  iigute  2  coiiesponds  to  the  mattix 
element, 

Af  (12;  34)  =<  nihitf-i:  njiao-j  U  I  nshso’s;  >4  .  (95) 

The  matrix  element  contains  four  terms  with  dilTetent  arrangements  of  the 
coordinates  and  the  wavefunctions.  For  rimplicity  of  notation  we  write  the 
siiQple  piodnct  states  as  1 12  >=  (ri)^n,k,(rj).  Then,  equation  (95) 

expands  as 

3f(12;34)=  (96) 

2  <  34  I  Vet  1 12  >  <  43  I  Vet  1 12  >  — 

<  34 1  Vet  1 21  >  <  43  i  1 21  >1  . 

(97) 

For  a  ^ven  initial  state  spin  there  are  three  distinct  phjrical  sitnations, 
with  equal  proability  of  occurence  corresponding  to  the  different  conAgura- 
tions  of  the  remaining  spiv  indices.  They  are: 

1.  ^1  S  O’]  =  0’s  = 

2.  Or  ^  O’]  :  =  ct,  O]  =  o’s 

3.  #  O’! :  =  «’4i  s=  o’t 

The  rate  for  each  configuration  must  be  calculated  separately,  and  then 
summed  to  give  the  total  rate  independent  of  spin.  iV>r  instance,  if  we 
define, 

Mx  =<34|K.|12>. 

Mi  =<43IV,«|12>, 

(98) 

then  the  squares  of  matrix  elements  corresponding  to  the  spin  configurations 
in  the  list  above  (the  probabilities)  can  be  written  in  terms  of  Afi  and  if] 

as: 

1.  \Mx-Mi\^ 

2.  lAfrl’ 

3.  |Af,|’. 

The  sum  of  these  probabilites  pves  the  square  of  the  total  effective  matrix 
element,  summed  over  all  internal  spins  for  a  given  initial  spin  which  we 
designate  as  Af^  write  as: 

M,\t  =  2\  Afi  I*  +  2|  Af]  I*  -  (Af;Af]  +  A/,*Afx).  (99) 


Thu,  we  need  only  to  coleolate  Mi  aad  i/j  to  find  the  toUl  rnte.  Further* 
more,  if  we  find  nn  expression  for  M},  Afj  is  ensilj  found  bj  exchanging  the 
final  state  indices,  so  we  onljr  calculate  Mi. 

The  simplest  way  to  include  frequency  dependent  screening  in  calculat¬ 
ing  Ml  above  and  the  associated  scattering  rate  5(13;  34)  is  to  calculate 
the  two-particle  propagator  from  |  12  >  to  j  34  >  in  time  t,  square  it  to 
obtain  the  probability,  and  take  the  time  derivative  for  long  times  (time 
rate  of  change  of  the  proability  of  going  from  1 13  >  to  1 34  >): 

5(12;  34)  =  ^lim  ^|<  34  |  K.  1 12  >|>  (100) 

This  is  wholly  equivalent  to  the  expression  in  section  2.  In  figure  2,  the 
electron  in  state  1  feels  the  screened  conlomh  potential  of  state  2  at  time 
t'  and  state  2  feels  state  1  at  time  The  dielectric  fnnction  retards  the 
effects  of  states  1  and  3  on  each  other;  for  iutanee,  even  if  dectron  2 
has  passed  electron  1,  the  valence  dectrou  may  still  ^  readjusting  to  its 
passage  and  this  can  effect  electron  1. 

We  write  Mi  —  <  34  |  )  12  >  by  letting  states  1  and  2  propagate 

&ecly  until  times  i*  and  t**,  and  then  scatter  into  states  3  and  4  which 
propogate  until  time  4.  The  potential  fSelt  at  each  time  is  weighted  by 
«~'(r''  —  —  t').  Then  we  mut  integrate  over  all  If  and  If.  This  is 

basically  a  restatement  of  the  Feynman  rules  for  a  first  order,  two-particle, 
time-dependent  interaction  which  |pve 

Mi^j^tj d*rVj(r)^;(rOl>r(r)^»(r')i; 

(101) 

with 


®  ®  (102) 
In  equation  (101)  we  have  written  the  Coulomb  potential  as  a  Fourier  series 
over  vectors  in  the  reciprocal  lattice  q.  If(t)  is  easQy  evaluated  by  making 
a  change  of  variables  in  the  t**  integral  to  t*  —  t".  Then, 


/,(t)  =  j  dtVt‘'*~^*J*‘*e’*(q,t''). 

Jq  Jf-t 

(103) 

Since  we  need  only  the  large  t  limit,  the  inner  integral  becomes  c~*(q,W}  — 
(•>4).  In  (101)  we  can  evaluate  the  spatial  coordinate  integrals  by  using  the 
trick  in  equation  (56)  section  4. 


=  ^53**«>**+*+®***»***-*+®* 


(104) 


In  the  above  equation,  Gi  and  (?>  ate  arbitrary  reciprocal  lattice  vectors. 
For  a  given  set  of  q  is  determined  up  to  a  tecipro^  lattice  vector.  We 
have  the  following  etyst«l  momentum  conservation  laws: 

1.  he  —  4j  +  q  +  Gi  s=  0 

2.  hs  —  hi  —  q  +  f7]  =  0 

3.  hi  +  hj  —  hs  •“  he  =  Go  =  Ci  +  Gj 

The  last  equation,  which  follows  from  the  first  two,  asserts  the  conservation 
of  crystal  momentum  for  the  entire  process.  Furthermore,  we  can  write  each 
of  the  Bloch  wavefunctions  as  sums  over  reciprocal  lattice  vectors  (since 
they  are  lattice  translation  invariant): 

t^»(r)  =  5;s»»(G)e*«'.  (105) 

a 

Then, 

Ml  =  5»,+*,.»,+»,+o,  (106) 

OtOtOt 

+  Ge  ~  Gi  —  Go)s|(Gs)aU^«) 

Vol  1  hi  —  hs  +  Gi  —  Gs  I* 

(107) 

This  form  is  particularly  well  suited  for  numerical  calculations  sbce  the 
expansions  in  (105)  can  be  calculated  by  psenoopotential  band  structure 
calculations.  To  calculate  the  scattering  rate  given  in  (100),  we  must  square 
(106)  and  tahe  the  time  derivative.  Tite  only  time  dependence  in  the  mag¬ 
nitude  squared  of  /,(<)  is  in  the  i*  integral  in  equation  (103).  The  time 
derivative  of  this  gives  the  Golden  rule  rate  expression  as  in  equation  (26). 

Now,  we  can  combine  (106),  (99)  and  (100)  to  obtain  the  total  impact 
ionisation  scattering  rate  5(12;  34)  summ^  over  secondary  particle  s^ns: 

5(12;  M)  =  ^  1  Mm  +  5,  -  J?,  -  .B4)6».+i..».+i.+o,(108) 

Mm  =  I  I’  +  I  h/,  I*  -  \(m;m2 + if; jf») 

Mi 

qi  —  hr  —  hj  +  Gi  —  Gs 
qj  =  hr  -  h^  +  Gt  —  G4 
Go  =  hr  +  hj  -  hj  —  h4 

This  expression  was  first  obtained  by  Kane  [38]. 

To  get  the  total  scattering  rate  we  need  only  snm  over  the  two  inde¬ 
pendent  h  vectors  and  the  secondary  particle  band  indices  in  in  5(12;  34). 
Given  nr  and  hr,  the  total  impact  ionisation  scattering  rate,  5«i,  from  that 


_  4-  <*4  **  —  <7o)^(f?s)s;(G4) 

0,^0,  V'olq?e(q,w,-W4) 


state  is  simplj, 

Su(niki)  =  53  535(12:34).  (109) 

The  wavevectoi  sums  (109)  run  over  the  first  BZ.  Go  is  then  the  Umklapp 
wavevector  necessary  to  ensure  that  hj  is  in  the  first  BZ  for  a  given  ht. 

To  gain  some  qualitative  insight  into  equation  (109)i  we  look  at  the 
two  conservation  laws  it  implies-conservation  of  crystal  momentum  and 
conservation  of  energy.  First,  energy  conservation  requires  the  initial  elec¬ 
tron  to  be  at  least  £f  (gap  energy)  above  the  bottom  of  the  conduction 
band  in  order  to  excite  a  valence  electron  across  the  gap  and  into  some 
conduction  band  state.  Therefore,  is  a  minimum  energy  for  impact  ion¬ 
isation.  Additionally,  the  conservation  of  crystal  momentum  farther  shifts 
this  minimum  energy.  The  smallest  energy  necessary  to  initiate  an  im¬ 
pact  ionisation  event  in  a  particular  band  structure  is  known  as  the  impact 
ionization  thrtzhold  energy  Eih‘ 

Some  authors  have  defined  a  wavevector  dependent  Btk  by  minimising 
Ei  (the  initial  electron  energy)  for  a  pven  kt  with  the  energy  and  momen¬ 
tum  conservation  constraints.  This  minimisation  procedure  requires  that 
(39] 

^hEzik)  =  Vfcf?s(Jfc)  =  VzE4(k).  (110) 

We  shall  re-evaluate  the  usefulness  of  this  result  later.  Now  we  will  look  at 
the  first  attempt  to  evaluate  the  impact  ionisation  scattering  rate  known  as 
the  Keldysh  formula.  The  Keldysh  formula  makes  two  rather  extreme  ap¬ 
proximations.  The  first  is  to  approximate  the  matrix  element  as  a  constant, 
and  the  second,  to  assume  simple  parabolic  bands.  Of  these,  the  second  is 
very  poor,  yet  the  resulting  formula,  known  as  the  Keldysh  formula,  has 
received  considerable  attention  in  the  past.  The  result  pven  by  Keldysh 
involving  the  two  adjustable  parameters  P  and  B  is  (40): 

SMEtk) 

Typically,  1  <  P  <  2  and  SfK  is  the  total  phonon  scattering  rate.  Although 
this  phenomenolopcal  form  is  simple  to  use,  the  approximations  leading  to 
it  disregard  the  important  features  of  the  band  structure,  which  at  the 
high  enerpes  involved  in  impact  ionisation,  has  no  resemblance  to  ample 
parabolic  bands. 

A  much  better  solution  for  the  impact  ionisation  rate  in  silicon  was 
given  by  E.  0.  Kane  [38]  in  1967.  He  numerically  calculated  the  rate  in 
equation  (109)  using  the  expressions  developed  in  equation  (108).  First  he 
calculated  the  pseudopotential  band  structure  and  wavefunctioiu  (see  the 
chapter  by  Fischetti  and  Higman)  and  used  them  to  calculate  the  dielectric 
function  in  the  RPA  approximation  (see  equation  (43)).  The  summations 
in  equation  (109)  were  evaluated  by  em^oying  a  Monte  Carlo  integration 
algorithm.  The  results  of  this  calculation  for  5m(E),  {Su{k)  averaged  over 
all  k  with  E(h)  =  E)  are  shown  in  figure  3. 

Figure  3  also  plots  the  impact  ionisation  scattering  rate  in  silicon  calcn- 
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Figoie  3:  Impact  iooifation  rates  in  s~*  for  silicon  nTetoged  over  initial 
electron  energy  measured  from  the  bottom  of  the  conduction  hand.  Dashed 
curve,  Kane's  result;  solid  curve,  no  collision  broadening;  dotted  curve, 
collifion  broadening. 
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Enasr  Measured  from  Cenduoinn-Biiid  (cV) 

Figure  4:  Secondaries  produced  by  impact  ionising  electrons  in  silicon. 
Solid  curve,  initial  electron  of  1.8  eV;  dashed  curve,  3.5eV;  dashed-dot, 
3.5eV;  dott^,  4.5  eV.  The  two  sets  of  carves  correspond  to  holes  for  the 
B  <  —Ea,  and  electrons,  for  £  >  0. 


Uted  to  include  the  eflecti  of  coUition  btoadening  and  the  btra-coUisional 
field  effect  (see  section  2)  [41].  The  effects  of  high  fields  and  high  phonon 
scattering  rates  on  the  impact  ionisation  process  shift  the  threshold  down 
fioin  that  obtained  by  Kane.  Figure  4  shows  the  distribution  of  secondaries 
(two  final  conduction  electrons  and  a  bole)  produced  initiating  electrons  of 
various  energies.  These  distributions  include  collision  broadening,  which 
has  a  large  effect  neat  threshold. 

Figure  5  addresses  the  question  of  h-space  anisotropy  for  the  impact 
ionisation  scattering  rate.  As  discussed  above,  it  is  possible  to  calculate  a 
wave-vector  dependent  threshold  which  would  seem  to  be  relevant  to  the 
anisotropy  of  the  scattering  rate.  It  is,  however,  difficult  to  guess  an  ap¬ 
propriate  form  for  the  scattering  rate  as  a  function  of  this  threshold  which 
is  consistent  with  the  true  physics  in  equation  (109).  For  instance,  the 
wave-vector  dependent  scattering  rate  may  not  even  be  directly  related  to 
these  thresholds  for  energies  greater  than  ^i\(h).  Furthermore,  if  collision 
broadening  is  included  (an  important  effect  for  high  energy  electrons  for 
which  the  phonon  scattering  rate  is  high)  the  threshold  condition  is  greatly 
relaxed.  Figure  6  shows  the  impact  ionisation  rate  in  silicon  for  electrons 
on  the  equi-energy  surfaces  E(Jt)  =  2.SeV  and  B(k)  —  S.OeV  for  electrons 
in  the  second  conduction  band  in  the  h«  =  0  plane  [41].  As  can  be  seen 
the  scattering  rate  shows  little  anisotropy.  In  addition  this  energy  range  is 
important  for  impact  ionization  in  tismsport  calculations  [42],  so  it  is  prob¬ 
ably  a  good  approximation  to  use  the  average,  energy  dependent  scattering 
rates  given  in  figure  3,  for  most  Monte  Carlo  simnlations. 
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Figure  5:  h-space  anisotropy  of  the  impact  ionisation  scattering  rate  in 
silicon  for  electrons  with  fcj  =  0  on  the  equi-energy  surfaces  E  =  2.5eV 
(solid  line)  and  E  =  S.OeV  (dotted  line).  The  angle  is  in  units  of  x  from 
the  hy  axis. 


6  Ionized  Impurity  Scattering 

Fot  even  modeiately  doped  lemiconductott,  ionised  imparity  scattering 
plays  an  important  role  and  can  dominate  the  total  scattering  rate  for  high 
doping  concentrations.  At  300K  when  doping  densities  teach  levels  above 
10^'cm*,  scattering  rates  for  ionised  donors  (acceptors)  become  comparable 
to  low  energy  phonon  scattering  rates,  while  at  low  temperatoies  ionised 
impurity  scattering  becomes  even  more  important  since  eqoilibrium  phonon 
populations  disappear  exponentially  as  T  -*  0.  Thus,  in  cases  of  high 
doping  or  low  temperatures,  an  accurate  Monte  Carlo  model  must  include 
ionised  imparity  scattering. 

Ionised  imparities  ate  usually  assumed  to  be  rimple  Coulomb  potentials 
with  a  charge  Zt  (e  the  electron's  charge).  IVpically  they  are  associated 
with  ionised  acceptors  or  donors.  At  first  glance  one  may  propose  that  the 
ionised  impurity  perturbation  K(r)  is  just  this  Coulomb  potential,  calcu¬ 
late  the  scattering  rate  for  one  impurity  and  multiply  by  the  total  number 
of  '‘^parities.  Unfortunately,  the  calculation  of  ionised  impurity  scattering 
rates  in  this  way  is  complicated  by  the  long  range  nature  of  the  Coulomb 
potential  associated  with  the  ions.  Unscreened,  the  interaction  of  a  single 
carrier  ud  an  ion  leads  to  a  diver^ng  scattering  rate  which  is  of  course 
an  unphysical  result.  Therefore,  the  view  of  the  ionised  impurity  interac¬ 
tion  as  one  electron  interacting  with  a  single  Conlombie  potential  cannot 
adequately  describe  the  situation.  Although  there  are  many  problems  with 
this  simplified  model,  the  two  most  fiagrant  are  the  neglect  of  mobile  charge 
screening,  and  the  correlation  of  the  other  charges  present.  Each  of  these 
limits  the  effective  length  of  the  Coulomb  interaction  and  removes  the  sin- 
gularity.  The  first  is  treated  by  the  Brooks-Herring  Model  [43],  and  the 
second  by  the  Conwell-Weisskopf  model  [44].  A  third  modd  proposed  by 
Ridley,  called.  Third  Body  Btclusion,  recon^es  these  two  approaches  [45]. 
We  treat  each  in  turn. 

0.1  The  Brooks-Herring  Model 

The  Brooks-Herring  model  assumes  that  the  electron  (or  hole)  interacts 
solely  with  one  ionised  impurity  site  and  deals  with  the  many-body  effects 
by  introducing  mobile  charge  screening.  If  mobile  charge  densities  are  high, 
it  u  necessary  to  include  the  effects  of  electronic  screening,  and  these  effects 
will  limit  the  scattering  rate  to  finite  values. 

Because  the  Brook-Herring  modd  assumes  that  a  carrier  sees  each  charge 
one  at  a  time,  in  a  semiconductor  sample  with  Np  ionised  donors  and  NJ 
ionised  acceptors  we  must  add  the  scattering  rate  due  to  each  separatdy 
instead  of  using  the  net  charge  present.  Thus,  compensated  ions  can  con¬ 
tribute  to  the  total  scattering  rate  as  well.  For  this  modd,  we  will  assume 
Ntot  total  ionised  impurities  of  charge  Ze. 

Next  we  turn  to  Thomas-Fermi  screening  theory  of  section  3,  equation 
(46).  If  the  distribution  function  is  Maxwellian,  (46)  reduces  to  the  Debye- 
Huckd  screening  formula,  introducing  a  screening  length  proportional  to 


the  iquaie  toot  of  the  mobUe  curiet  deniity  n: 


= 


*0€,K,T  * 


The  Fourier  tiansform  of  the  scteened  loniied  impuritjr  potential  1^(9)  is: 


V,(«)  = 


e,eofl’  «r^(«) 


Ze 


For  Kfr  0  in  equation  (113)  thete  is  no  diveigence  fot  small  9,  and 
theiefore,  there  is  no  problem  with  inUnite  scattering  rates.  For  degener¬ 
ate!/  doped  semiconductors,  the  situation  is  more  complicated,  however  an 
effective  screening  length  can  still  be  calculated  from  (48). 

In  teal  space,  the  Coulomb  potential  acquires  an  exponential  tail  as: 


Vi(r)  = 


4xe,eo }  r  -  II  | 


where  R  is  the  location  of  the  ion.  It  should  be  emphasised  that  the  screened 
potential  approrimation  is  not  necessary  to  limit  the  cross  section  to  finite 
values.  In  fkct,  at  low  mobile  charge  density,  kfr  may  not  be  sufficient  to 
screen  the  potential  and  a  model  similar  to  the  Conwell-Weisslcopf  method 
must  be  used. 

In  this  section  we  will  assume  rimple  plane  wave  states  for  the  riectrons 
which  is  a  good  approximation  fot  low  energies.  Ionised  impurity  scattering 
is  essentially  important  only  for  low  energy  electrons,  because,  the  phonon 
scattering  rate  dominates  the  impurity  scattering  rate  fot  high  diectron 
enetpes. 

We  write  the  perturbation  matrix  element  between  plane  wave  states 
as  Mkv  where  i  is  the  incident  wave  vector  and  h'  is  the  scattered  wave 
vector  as 


=  7/*-“' 


4xe,<o  I  r  -  H I 


(Vo(,  the  volume  of  the  crystal,  enters  because  of  normalisation).  The 
matrix  element  in  equation  (115)  is  found  by  talcing  the  integral  over  into 
spherical  coordinates  with  the  asimnthal  axis  being  in  the  direction  of  9  = 
y  —  k.  The  result  is: 


Mik<  — 


iZt*  *  9 

e,co«  Vd  -f-  9* 


Fbally,  the  scattering  rate  5(h,  h')  is  pven  by  the  Fermi  Golden  Rule  (see 
equation  (38)).  Fkom  equation  (116)  we  can  calculate  the  total  scatter¬ 
ing  rate  Stat  by  integrating  S{Jk,V)  over  all  V.  Assuming  spherical  non- 


puabolie  bands, 


Sttt  = 


Vol 


2x  A* 

li  (1  +  2af?/)m«  8r*  " 

r  r  k^dmE;-E,) 
J~i  Jo  (i>  +  *'»  -  24'*«  + 


(117) 


Here,  ti  is  the  cosine  of  the  angle  between  k  and  k’.  Since  the  ion  is  strongly 
coupled  to  the  lattice  and  is  much  mote  massive  than  the  electron,  the 
interaction  can  be  treated  as  elastic,  and  we  can  take  the  term  Ej  equal  to 
the  energy  of  the  final  electronic  state  with  wave  vector  V  and  Ei  equal  to 
the  energy  of  the  initial  electronic  state  with  wave  vector  k. 

The  total  scattering  rate  is  [43] 


Stai 


^  _ Ni_ 

32x(2m’ 
(1  +  2aE) 
^l\E) 


l{kTr/2k)m  +  (kr^/2ife)»)) 


(118) 


where  ffj  is  the  total  density  of  impurities  {Nx  —  IftotlVol)-  It  is  important 
to  notice  that  we  have  simply  summed  the  interaction  of  a  plane  wave  and 
a  Coulombie  center  over  all  the  centers  present;  we  have  not  calculated 
any  many*body  effect  other  than  mobile  charge  screening-the  calculation 
is  a  two-body  calculation.  The  probability  of  a  scattering  into  angular 
increment  d$  about  $,  (0  is  the  angle  between  k  and  k*)  P{0)dS  is 


P{0)dB  s 


nn9d$ 

(2*>(1  -  cos«)-b*»,)‘ 


(119) 


To  relate  P[0)d0  to  a  uniformly  chosen  random  number,  we  need  only  find 
the  normalised  probability  that  9  lies  between  0  and  0.  The  result  of  this 
normalisation  aUows  ns  to  determine  9% 

^  -  1  +  4r(*/*«.)>  • 

Thus,  we  can  stochastically  determine  the  final  state  with  only  two  random 
numbers,  one  for  ^  and  one  for  9.  kf  must  equal  k  because  the  scattering  u 
elastic  and  k  is  predetermined.  If  the  bands  are  ellipsoidal,  non-parabolic 
bands,  we  can  use  the  Herring- Vogt  transformation  mentioned  in  section  4 
equation  (84).  We  need  only  make  two  modifications.  First,  we  replace  m* 
by  mj)  in  equation  (118).  To  choose  the  final  state  we  pick  a  k'  using  (120). 
Then,  we  transform  it  back  &om  the  Herring- Vogt  space  (see  equation  (84)). 

It  should  however,  be  emphasised  that  the  screened  potential  is  not 
necessary  to  limit  the  scattering  rate  to  finite  values.  In  (act,  for  low  mobile 
charge  density,  kn  may  not  be  sufllcicnt  to  screen  the  potential  and  a 
model  such  as  the  Conwell- Weisskopf  model  may  be  more  correct. 


6.2  Conwell-Weisskopf  Model 

For  comparison  with  the  Biooks-Herring  model,  we  compute  the  scattering 
cross  section  for  the  Conwell-Wcisskopf  model  which  assumes  screening  by 
other  ions  instead  of  screening  by  mobile  charges.  The  Conwell-Weisskopf 
model  assumes  that  the  electron  is  (airly  well  localised  instead  of  being 
an  infinite  plane  wave  as  assumed  in  the  Brooks-Herring  model.  This  as¬ 
sumption  remains  useful  up  to  moderate  energies  where  the  electron  is 
considered  a  wave  packet  following  a  classically  defined  orbit.  By  localising 
the  electron,  we  can  see  that  in  certain  positions  (midway  between  ions), 
the  Coulomb  forces  from  each  tend  to  cancel  and  there  is  no  interaction. 
In  this  spirit,  Conwell  and  Weisskopf  limit  the  impact  parameter  b  to  half 
of  the  mean  distance  between  ions, 

6m.,  =  (121) 

where  the  impact  parameter  for  classical  Rutherford  scattering  is  the  closest 
approach  of  the  electron  to  the  ion  if  its  path  weren’t  deflected  by  the 
Coulombic  force.  Thus,  this  model  effectively  screens  out  the  Coulomb 
force  to  the  average  distance  between  ions,  and  the  unscreened  Coulomb 
potential  can  be  used  to  calculate  the  scattering  rate  classically  with  the 
Rutherford  scattering  model  [44].  The  scattering  rate  evaluated  in  this  way 

5..i(ifc)  =  (122) 

where  o(k)  is  the  electron’s  velocity.  The  electron  equation  of  motion  states 
(hat  r(A)  =  (l/h)VkE{k).  Pot  parabolic  bands  then,  v{k)  =  (hJt)/m*. 

Not  only  is  the  Conwell-Weisskopf  model  more  applicable  for  low  doping 
densities  than  the  Brooks-Eerring  model,  but  it  also  behaves  more  reason¬ 
ably  for  high  doping  densities  than  the  Brooks-Eerring  model.  E  we  plot 
mobility  versus  Nj  for  both  models,  we  see  that  the  two  agree  very  well  for 
Nj  <  10^^  cm*. 

6.3  Third  Body  Exculsion 

Both  models  presented  above  imply  two-body,  nearest-scatterer  processes, 
but  don’t  expressly  prohibit  scattering  &om  more  distant  scattering  centers. 
Ridley  argues  that,  for  the  sake  of  consistency,  the  differential  cross-section 
ff{k,  8, 4)  must  be  weighted  by  the  probability  the  the  scattering  is  a  nearest- 
scatterer  process.  In  doing  so,  he  has  found  that  the  Umiting  cases  of  both 
the  Brooks-Eerring  and  the  Conwell-Weisskopf  models  can  be  obtained. 
The  method  of  weighting  the  cross-section  in  this  manner  is  called  third- 
body  exclusion.  (Note:  the  classical  scattering  cross-section  is  related  to 
the  scattering  rate  by  5(k,^,0)  =  Nj  v{k)ff{k,8,^)). 

In  order  to  determine  the  appropriate  factor  we  again  employ  the  classi¬ 
cal  notion  of  the  impact  parameter  6.  Ridley  has  calculated  the  probabUity, 
P{b),  that  no  scattering  center  exists  with  impact  parameter  less  than  h 


from  the  propability  p  =  ixNiah  ib  that  inch  a  center  exists: 

P(6)  =  .  (123) 

Here,  a  U  the  average  distance  between  ions.  Therefore,  to  prohibit  third- 
body  processes,  the  probability  of  a  scattering  event  occurring  in  the  solid 
angle  d(l,  mast  be  mnltiplied  by  P(b).  Thns,  we  need  only  calculate  h 
corresponding  to  the  cross-section  o  we  wish  to  nse.  If  we  take  o  from  the 
Brooks-Herring  model,  we  can  obtain  the  limiting  cases  of  both  models,  so 
we  follow  Ridley  and  calculate  the  corrected  diflerential  cross-section  for  it 
(45). 

9c  the  corrected  cross-section,  can  be  calculated  easily  from  the  Brooks- 
Herring  cross-section  caa  hy  the  following: 

ec  -  2x  J’ ^6de  (124) 

Roer  and  Widdershoven  give  the  scattering  rate  Sc  —  Jfjv(k)9c  ss  [46] 


This  scattering  rate  incorporates  screening  of  both  types  (mobQe  and  fixed 
charges)  because  it  contains  both  kxa  ^  ionic  screening  cutofT  from  the 

term  in  equation  (123).  Furthermore,  it  has  the  advantage  that  the  total 
scattering  rates  are  ten  to  a  hundred  times  lower  than  the  peak  valnes  of 
the  Brooks-Herring  and  Conwell-Weisskopf  models.  This  is  quite  important 
for  Monte  Carlo  nmulations  because  the  higher  the  scattering  rate  for  all 
mechanisms,  the  smaller  the  time  step  must  be,  forcing  simulation  runs  to 
take  much  longer  times.  Roer  and  l^ddershoven  also  show  that  at  least  in 
the  case  of  G^ks,  the  low  field  imparity  limited  mobilities  agree  weO  with 
experimental  results  [46]. 

To  find  an  expression  for  the  angular  distribution  from  the  third-body 
exclusion  modd,  we  notice  that  the  probability  that  an  electron  has  an 
impact  parameter  h  is  pP(b)  —  db.  This  is  simply  the 

pr^uet  of  the  probabOity  that  a  scattering  center  exists  at  a  distance  b 
and  the  probability  that  no  other  scattering  center  is  closet.  We  can  relate 
P(b)db  to  a  uniform  random  number  r  :  r  €  [0, 1]  by  the  relation: 

r  =  .JsLll - =  —s - -T  •  (126) 

jJ— -  1 

Since  the  relationship  between  r  and  b  is  transcendental,  the  Von-Neumann 
rejection  method  must  be  used  to  select  a  value  of  6.  Then  the  scattering 
angle  8  can  be  determined  from  b  by  inverting  the  relation 


Th*(d)  =  J’  9BB(kt8')sin9'<U' 


(127) 


Then,  the  lelntion  between  6  and  b  U 


1 

/2k*)r6» 

.  1  ] 

2k> 

\ 

(128) 


with 


,  ^  (Ze*yil  +  2aE)m-i 


(129) 


Although  (till  n  rather  crnde  model,  third-bodjr  ezclnsion  is  a  com* 
promise  inclnding  screening  of  both  types  which  limits  scattering  rates  to 
manageable  levels  and  fits  experimentally  obtained  low  field  mobilities  fairly 
well.  At  high  fields,  the  Born  approximation  breaks  down  and  other  meth* 
ods  mast  be  sought  to  include  collision  broadening  effects. 
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Appendix 

Table  of  Physical  Constants 
and  Phonon  Scattering  Parameters  for  Silicon 

The  following  table  gives  typical  parameter  values  and  material  con* 
stants  for  Si  which  occur  in  the  fotmnlas  of  section  4.  Note  that  intervalley 
phonon  energies  are  given  in  degrees  Kelvin  (h(i»  s  hj|7^a«nM)-  Also,  the 
models  shown  here  are  two-valley  models-  X  and  L  valleys. 


1  Units 

IKSaiSI 

iiomiaji 

IK23CSI 

a 

X 

S.43 

5.43 

5.43 

P 

g/cm* 

2.329 

2.329 

2.329 

2.329 

« 

<0 

11.7 

11.7 

V, 

10*  cm/s 

9.04 

9.037 

9.04 

9.04 

tn.1 

mo 

0.19 

0.1905 

0.19 

0.19 

m,t 

mo 

0.9163 

0.916 

0.916 

etm 

•  eV 

0.5 

0.5 

Xi, 

eV 

9.5 

KM 

9.0 

9.5 

mil 

mo 

0.12 

0.12 

0.12068 

jTta 

mo 

1.59 

IBBI 

1.59 

1.5942 

1  X-X  Intervalley  Scattering  | 

eM 

220(0.3) 

19 

550  (1.9) 

685(2.0)  1 

630  (4.0) 

630  (4.0) 

685  (1.9) 

140  (0.5) 

215(0.8} 

210  (0.8) 

215  (1.1) 

ttiltUUll 

700(3.0)  1 

EIIIQDII 

720  (4.3) 

1 _ ^ 

'j  Intervalley  Scattering  | 

fm 

1  :€ 

672  (4.0) 

672  (1.8) 

■9Mb 

634(2.0) 

634  (4.0) 

634  (1.8) 

19 

liiSl 

480  (1.8) 

480  (1.8) 

■■Bi 

197  (1.8) 

197  (1.8) 
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